Anomalous Hall effect from frustration-tuned scalar chirality distribution in Pr2lr207 
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We analyse the Ising Kondo lattice model on a pyrochlore structure in order to study the anoma- 
lous Hall effect due to non-coplanar magnetism. We focus on the frustration-induced spatial inhomo- 
geneity of different magnetic low-temperature regimes, between which one can efficiently tune using 
an external magnetic field. We incorporate non-magnetic scattering on a phenomenological level so 
that we can distinguish between the effects of short-range correlations and short-range coherence. 
We obtain a Hall conductivity (fn) as function of field strength and direction which compares well 
to the experimental data of Pr2lr207. In particular, we show that the observed peak in cth for 
H II [111] signals the crossover from zero-field spin ice to Kagome ice. 

PACS numbers: 71.10.Fd, 71.23.-k, 71.27.+a, 72.10.-d 



The properties of itinerant degrees of freedom on geo- 
metrically frustrated lattices are only poorly understood. 
One promising avenue for studying the interplay of frus- 
tration and itinerance are hybrid systems where itinerant 
electrons interact with localized magnetic moments sub- 
ject to strong frustration. The latter can exhibit various 
exotic phases incorporating peculiar spatial correlations 
[l[ . It is therefore natural to ask whether these bequeath 
their unusual behavior to the itinerant electrons, result- 
ing in novel types of behavior for the composite system. 

The anomalous Hall effect (AHE) is one particularly 
striking resulting phenomenon AHE was originally 
associated with ferromagnetic conductors with strong 
spin-orbit interaction [3|-l5|- However, AHE has recently 
been reinterpreted in a broader context, including non- 
coplanar magnets as promising candidates for its emer- 
gence [1, H-li]. 

Prominently, the compound Pr2lr207 shows a unique 
Hall response. It is composed of two interpenetrating 
pyrochlore lattices. Ir 5d electrons form a conduction 
band on one, while the localized Pr 4/ moments reside 
on the other, and develop spin-icc-type correlation at low 
temperature[l, It is quite plausible that the spin 

scalar chirality of the spin ice manifold gives rise to non- 
trivial features in the Hall response, particularly strik- 
ingly in zero field [l2j • In addition, the observed Hall con- 
ductivity is highly anisotropic and non-monotonic, with 
a prominent peak around H ^ 0.7 Tesla for H || [111] 
B 13 . [l^ . P ioneering analyses of the pyrochlore conduc- 



tors |14l . Il5| have considered spatially periodic structures 
for the localized moments. It is now natural to ask how 
spatial aperiodicity - arising from the geometrical frus- 
tration of the spin-ice local moments - manifests itself in 
the nontrivial Hall response observed in this compound. 

We analyse the dependence of a^y on the coupling be- 
tween itinerant and local degrees of freedom; field di- 
rection and strength; as well as damping which we in- 
clude for the itinerant electrons to phenomenologically 
take into account non-magnetic scatterers. Our results 



compare well with the experimental data for Pr2lr207 
with H II [100] and [111]. In particular, we find that 
the prominent peak observed for H || [111] can be at- 
tributed to the crossover from the zero- field spin ice state 
to Kagome-icc state: the latter is a state with perfectly 
field-aligned spins on the triangular layer with the other 
spins disordered but subject to the ice rule constraint [Til] . 
see Supplementary I. 

In the remainder of this work, we first present the 
model and analysis method for the relatively simple 
Kagome ice model. This is then extended to include 
the features necessary to make detailed contact with the 
Pr2lr207 data. 

The Kagome ice model describes localised non- 
dynamical Ising moments on the Kagome lattice [Fig. 
[1] (a), (b)] interacting with itinerant electrons through 
local fields at each site: The Hamiltonian is given by 
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The sum is taken over the nearest-neighbor (n.n.) 

sites. We simply choose hi = JSi, with exchange cou- 
pling J. The localized spins {Si} are subject to local 
easy-axis anisotropy, i.e. Si = ryiDi [rji — ±1), with 
D, = ^[1, -1, 1], ;^[1, 1, -1] and ^[-1, 1, 1], if I be- 
longs to sublatticc A, B and C. Here, we take a quenched 
average in terms of {Si} by imposing "Kagome ice rule", 
namely we impose for each triangle X^ieA = 

This allows a macroscopic number of spin configura- 
tions [l7, 18 1, including a uniform configuration [Fig. [T] 
(a)], and a huge number of disordered configurations [Fig. 
[l](b)]. Crucially, for each and all of these, the spin scalar 
chirality is uniform Sa • (Sf, x Sc) = Kq = — 4/3-\/3 for all 
the upward and downward triangles We can thus 

examine the effect of spatial disorder on AHE, while pre- 
serving uniform spin scalar chirality. 

For the calculation of Hall conductivity ctj-j,, we ran- 
domly generate a series of spin configurations under the 
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Hereafter, we set t 



7/1 = 1. 



FIG. 1. (color online), (a) A uniform configuration and (b) a 
representative of the disordered configurations satisfying the 
Kagome ice rule. Blue (red) arrows show spins Si correspond- 
ing to rji = 1(— 1). An example of the graph belonging to 
G[7] (see the main text) is shown in (b) with combined three 
dashed arrows, (c) Structure of double pyrochlore lattice. 
Green (Gold) tetrahedra constitute the Ir (Pr) pyrochlore lat- 
tice. Sublattice indices A, B, C and D are shown for the Pr 
lattice, for which an example of a spin ice configuration is 
shown, (d) One Ir tetrahedron surrounded by 4 Pr tetrahe- 
dra. Each Ir ion (i) has 6 neighboring Pr ions (jn ■ ■ ■ jis) 
forming a hexagon, as highlighted by a thick line. 



Kagome ice rule, {s|''''}. For each {S^^^}, the Hall con- 
ductivity is given as 



m.m' 

liii{{m.\Jx\m'){m'\Jy\m)) 

(2) 

by Kubo formula. Here, |m) and are the eigenenergy 
and corresponding eigenstate of Hamiltonian eq. ([l}. 
f(E) is the Fermi distribution function at zero temper- 
ature. Jx{y) is the x{y) component of the current op- 
erator, and V is the total volume of the system. Here, 
we introduce the phenomenological damping rate 1/r to 
take account of the finite lifetime of electrons due to non- 
magnetic impurities. While the magnetic disorder itself 
causes damping, non-magnetic scattering plays another 
important role in Hall conductivity, l/r sets a coherence 
length of electrons, which determines the effective spa- 
tial scale of spin scalar chirality. The Hall conductivity 
Gxy can be obtained after taking the configurational aver- 
^' ~ri;({S-^^}). We typically choose 



N. 



In Fig. [2ja) and (b), we show the dependence of axy 
on particle density, n, in the uniform [Fig. [TJa)] and 
disordered configurations [Fig. (Ub)], at relatively large 
damping, l/r = 1.0. In the uniform case, Gxy can be ob- 
tained as a summation of Berry curvature in momentum 
space, showing a steep change around n ^ 1/3, where the 
Dirac cone in the non-interacting band gives a singular- 
ity in Berry curvature. This singular behavior is absent 
in axy in the disordered case. 

J and l/r dependence of axy is summarized in Fig. 
^c) at n = 0.0977. Firstly, for small J, a cubic law, 
axy oc J^, is found in both uniform and disordered cases. 
This cubic law can perturbatively be ascribed to the mul- 
tiple scattering from triplets of spins exhibiting finite 
scalar chirality [23|- With increasing J, deviation from 
the cubic law is found at J ~ l/r. In particular, in the 
uniform case, axy becomes insensitive to l/r, and an- 
other scaling law, axy oc J, appears, suggesting the axy 
is described in terms of the Berry curvature in this region 
For J ^ t, the system falls into a double-exchange 
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regime: the itinerant electron spins are aligned with the 
localized spins, and axy takes on values only weakly de- 
pendent on l/r. 

In general, axy takes considerably different values be- 
tween the disordered case (a^y) and the uniform ordered 
one, a^y, as shown in Fig. ^d), where we plot the l/r 
dependence of axy at J = 0.05. This is most pronounced 
for small damping, l/r ^ J, where cr"y saturates, but 
the difference persists all the way to l/r 3> J. 

A pcrturbative treatment in h sheds light on the origin 
of difference between the two cases. To third order 123, 



^xy 



(41,42, 



(hi, X h,^)Wxyiii,i2,i3), (3) 



age, as a^ 

TVs = 100, and system size iV = 32 x 32 x 3 = 3072 sites. 



where the summation is taken over the iV(A^— 1)(A^— 2)/6 
triplets of sites (11,12,13), see Supplementary material 
II. This gives the Hall conductivity as summation over 
the triplets' spin scalar chirality with weighting factor 
Wxy{ii,i2,h)- It is instructive to resolve the Hall con- 
ductivity in the form of a graphical series expansion 

as axy = X] '^xy^ ■ Here, cri™^ is the total contribution 

m— 3 

from the triplets (11,12, is) belonging to the set of graphs 
G[m] composed of three segments of total length m (An 
example of a triplet € G[7] is shown in Fig.[TJb)). 

Fig. [5Je) and (f) show the partial summation S*!™' = 
J^iLs '^xy and the averaged weighting factor at each m, 
Wxy \ at J = 0.1 and l/r = 0.5. In the disordered 
case, different contributions for graphs of a given to > 5 
come with an effectively random sign, hence canceling 
against one another. By contrast, for the uniform case, 
the summation continues to oscillate until damping de- 
stroys coherence at larger to. This illustrates the different 
roles played by loss of correlations of the local moments 
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ant electrons (cia) on the Ir sublattice, and localized 
Pr moments ({S[j]}) on the other. The localized mo- 
ments are subject to Ising anisotropy as S[j] = VjlDj, 
with = -i=[l,l,l],-i=[l,-l,-l],-i=[-l,l,-l] and 

-^[—1, —1, 1], a j belongs to sublattice A, B, C and D, 
respectively. Each site i on the Ir sublattice has 6 neigh- 
bors (jii,Ji2 • • • , Jie) located on the honeycomb ring of 
the Pr sublattice, as shown in Fig. [1] (d). To describe 
1 -0 the interaction, we adopt the same Ising-Kondo-lattice- 
type Hamiltonian eq. ([T|) with the local field the sum of 
the 6 neighboring localized moments = ^[ji']- 
While the spin configuration can be determined from 
the equilibrium condition of Hamiltonian ([1]), we start 
from the phenomenology that the spin part is described 
by the nearest-neighbor spin ice Hamiltonian (j4|), in order 
to examine how spin ice correlations in Pr2lr207 affect 
its Hall conductivity. 
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FIG. 2. (color online). Dependence of axy on particle density, 
71, for (a) the uniform configuration and (b) the disordered 
configuration at 1/r = 1.00. (c) J dependence of axy at 
n = 0.0977 for the disordered case (dots) and the uniform 
case (solid lines), (d) The 1/r dependence of a^y at n — 
0.0977 and J — 0.05. The dashed line is a guide to eye. 
(e) Partial summation of Hall conductivity Si^K (f) The 
weighting factor VFi™' averaged over the graphs at each m. 



and loss of coherence of the itinerant electrons in the two 
respective cases. 

Now, let us apply this 3rd-order perturbative scheme 
to Pr2lr207. Since the exchange couphng J stems 
from the superexchange process between Pr and Ir ions, 
it is reasonable to assume J/t <C 1. We consider 
a double-pyrochlore lattice: two interpenetrating py- 
rochlorc lattices, as shown in Fig. [T] (c), with itincr- 



J. 



spin VjVj' 



H 



We use ?^spin in a standard equilibrium Monte Carlo sam- 
pling to obtain = 100 sets of {S[j]}. Though we are 
interested in the region T — 0, we introduce temperature 
T as a phenomcnological parameter to mimic the devia- 
tion from ideal spin ice due to the long-range RKKY in- 
teraction in the actual compound, and set T/Jspm = 0.5. 
Hereafter, we focus on the field directions H || [100] and 
[111] [2^. We set magnetic coordinates || [010] and 
By II [001] for H II [100], and e^, || [110] and ey || [112] for 
H II [111], and calculate cth = (Jxy We consider the low 
density region 0], and fix the particle density at n = 0.01. 
As system size, we adopt = 12 x 12 x 12 x 4 = 6912 
sites. 

In Figs. O (a) and (b), wc plot the magnetic field de- 
pendence of (Th at 1/r = 5.0 and 0.5. For extremely 
large damping, 1/r = 5.0, only the smallest triangles 
contribute to cth [Fig. WiA)] ■ In this region, the sign of 
(Th becomes opposite between H || [100] and H || [111], as 
expected in Ref. Q on the assumption of the local limit. 
However, the full magnetic field dependence of cth, espe- 
cially the low-field negative linear response in this local 
limit considerably deviates from the experimental results 
01 [Fig. El^a), inset]. The negative linear response comes 
from the large negative contribution at m = 3 [Fig. [3] 
(d)]: solely short-range spin ice correlations within the 
four-Pr cluster [Fig. [l](d)] do not give correct cth- 

In contrast, for intermediate damping, 1/r = 0.5, (Th 
shows quite similar behavior to the experimental data 
[3]. For small H, cth shows positive linear response irre- 
spective of the field direction [Fig. IS^b), inset]. In Fig. 
[3] (e), we plot the graph-resolved Hall conductivity at 
H = 0.4 for H II [111]. This plot shows that the spatially 
extended scalar chirality beyond the local limit {m < 10, 
or roughly 35A) plays a crucial role in the positive linear 
response. 
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FIG. 3. (color online). The magnetic field (H) dependence 
of an at n = 0.01 and J = 0.1 for H || [100] and [111] at 
(a) 1/r = 5.0 and (b) 0.5. The results for H || [100] ([111]) 
are plotted with open squares (diamonds). The insets are 
enlarged plots around H = 0. (c) H dependence of the prob- 
abilities P22, P31, and Pa- The crossing of P22 and P31 is 
shown in the inset. The vertical dashed line in (b) and (c) 
is a guide to eye. (d) (e) Graph-resolved Hall conductivity 
a^"^ and its partial summation S^-' at H/Jspin ~ 0.4 and 
H II [111] for (d) 1/t = 5.0 and (e) 0.5. (f) an of the muhi- 
orbital model. 



The low- field peak for H || [111] is tlie most conspicu- 
ous feature of the Hall conductivity in Pr2lr207. In pre- 
vious studies 0, [1^ , this peak is attributed to the spin 
flip crossover from the low-field spin ice state with dom- 
inant 2-in 2-out configuration to the high-field saturated 
state with 3-in 1-out and 1-in 3-out spin configurations. 
However, our analysis suggests a different picture. In 
Fig. [3] (c), we plot the probabilities that each tetrahe- 
dron is occupied by 2-in 2-out configuration (P22), and 
3-in 1-out or 1-in 3-out configurations (P31). This plot 
shows that P22 and P31 are almost constant at low fields, 
until the spin flip crossover happens at much higher field 
ff-6J,pi„ [H [Fig. He), inset]. 

The peak of cth seems rather related to the crossover 
from the zero-field spin ice state to Kagome-ice state. 
In Fig. [3] (c), we plot the probability that a spin on the 
sublattice A aligns parallel to the field (Pa), as an in- 
dicator of the Kagomc ice state. Pa changes from 0.5 
at 7? = to 1.0 at the Kagomc ice state. The peak of 



cth corresponds to Pa ~ 0.75, i.e. the midpoint of this 
saturation process, clearly showing that the peak signals 
the crossover to a Kagome ice state. Within the nearest- 
neighbor spin ice model used here, the crossover occurs 
at 7? ~ T. Accordingly, the peak is located at 7J ~ O.ST, 
see Supplementary material I. Indeed, the [111] magneti- 



zation takes M = M„ 



has a peak, with Mgat 
netization [1]. This Mp 
midpoint of saturation process of Ma 



^Mgat experimentally, when cth 
' 1.5/iB/Pi', the saturated mag- 
coincides with the value at the 
and smaller than 



TMsat expected at the spin flip crossover. 



Here, let us turn to the physical origin of the peak. 
As the magnetic field is applied, net spin scalar chiral- 
ity, and hence cth, is enhanced. On the other hand, the 
evolution to the Kagome ice state can be regarded as a 
partial ordering process of sublattice A, so that spatial 
disorder is reduced, suppressing the interference between 
the graphs resulting in a suppression of cth, as discussed 
above for Kagome ice model. It is tempting to note that 
the balance between the two effects gives rise to a promi- 
nent peak during the evolution to the Kagome ice state. 
Indeed, the peak becomes more prominent as 1/t is fur- 
ther reduced (not shown), reinforcing the subtle balance 
between the two competing effects. The sensitivity to 
1 /r may be confirmed from the systematic study of sam- 
ple dependence of (Th and resistivity. Further analyses 
are clearly desirable to elucidate this point. 

We finally comment on the quantitative aspect of our 
theory. Although our result is based on a single or- 
bital tight-binding model, the overall features of cth are 
insensitive to band structure. However, the amplitude 
of the Hall conductivity is sensitive to 'details' of the 
band structure. By adopting a realistic band struc- 
ture based on a multi-orbital tight-binding model with 
three t2g orbitals, we could obtain the Hall conductivity 
~ 30ri^^cm^^ at high fields, comparable to experimen- 
tal data [Fig. 131(f)] j with reproducing various features of 
experimental results; See supplementary material HI. 
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Supplementary Items for Anomalous Hall effect from frustration-tuned scalar 

chirality distribution in Pr2lr20'j 

I. CROSSOVER BETWEEN THE SPIN ICE AND THE KAGOME ICE STATE 




FIG. 1. (color online). A schematic picture of the crossover from zero-field spin ice to Kagome ice. Only the Pr lattice is 
shown. In Kagome ice the spins are field-aligned on sublattice A indicated by shaded ellipses, which form layered triangular 
lattices. 

The spin ice state is characterized by the so-called ice rule: out of the four spins on each tetrahedra, two point 
inward while the other two point outward, satisfying the 2-in 2-out constraint [Fig. [T]. Under the magnetic field 
H II [111], the spin ice state changes to the kagome ice state, which has perfectly field- aligned spins on the triangular 
layer with the other spins disordered but subject to the ice rule [Fig. [T]. Looking from the [111] direction, a pyrochlore 
lattice can be considered as the alternate stacking of kagome and triangular lattices. Since the spins on the triangular 
layers have easy axis parallel to [111], these spins tend to align prior to those on the kagome lattice under the magnetic 
field H II [111]. 

In Fig. [21 we plot the magnetic field (H) dependence of magnetization (M) of the nearest-neighbor spin ice model 
eq. (4) in the main text, obtained at T = 0.5Jspin- For H <T, M grows almost linearly, and reaches M = |Msat at 
H ^ T, with Mgat the saturated value of magnetization for H || [111]. For T < H < GJgpin, the kagome ice state is 
stabilized with nearly constant magnetization M ~ |Msat. For H > 6Jspin, the spins on the kagome layers also align 
in the field direction. As a result, the magnetization gradually approaches its saturated value, Al Mgat- 
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FIG. 2. (color online). The magnetization process of the nearest-neighbor spin ice model eq. (4) in the main text, obtained at 
T = O.SJspin. The magnetization (M) shows different behaviors at each magnetic field (H) regions. 



II. THIRD-ORDER PERTURBATION THEORY OF ANOMALOUS HALL CONDUCTIVITY 

In this section, we give a detailed derivation of the pcrturbative formula of Hall conductivity [eq. (3) in the main 
text]. We start with the Ising Kondo Hamiltonian [eq. (1) in the main text], in momentum space representation. 



k,s a, a/ 



(5) 



k,k' ss' a 



Here, we have introduced sublattice indices a and a'. is the number of unit cells. T-Laa'O^) is the kinetic energy in 
momentum representation. For example, for the Kagome ice model, it is given by 



— 2 cos(k • (ga — bb)) — 2cos(k • (ga — ec)) 

— 2 cos(k • (gb — ga)) — 2cos(k • (gb — Gc)) 

— 2 cos(k • (gc — ga)) — 2cos(k • (gc — Gb)) 



(6) 



Here, ga, gb and gc are the internal coordinates in a unit cell, corresponding to sublattice A, B and C. We consider 
the current correlation function 



(7) 



defined with the current operator 



_ dUgg' (k) I 



k,a,a' .s 



(8) 



and Matsubara frequency LOq ~ 2TTqT [q: integer). The Hall conductivity a^'^ can be obtained as an asymmetric part 
of the derivative of correlation function as 



d g^-(c. + ^0)-Q-A.(^ + ^o) 1 

^ = ^(<3 ' (27riT) - {2mT)). 



idu! 



(9) 
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We expand Q^'^{iujq) to third order of "H', 

l/T 1/T 1/T 1/T 

'^'^ f 1 iu),.T f 1 f 1 f 1 dHab dHa'b 



, f , f , \ ^ \ ^ dHab dHa'b' \ ^ \ ^ 1 \ ^ 1 \ ^ 1 \ ^ 



dr e*'^'"' 

6V J J "J '"J ^ , dkt" dk" ^ ^ N ^ N 

Q ksab k's'a'b' a^a^a^ ap^ kik^ ^2^2 ^3*^3 

w n hi hi cr" <J^ ,a'' , 

Kl— K^-Qi k2— k2,a2 ks— k3,a3 SiS-,^ -5252 ^^^3 

X ('^4asMckt.,(r)4,^,,,(0)ck'h's'(0)4^^^,^(rl)c^alsi(n)42a2^2(^2)c^a2S^(T•2)43a3S3 

4 J \ ^ \ ^ \ ^ dT-Lab dHa'b' \ ^ 1 u \ 

= 'VIV^ -Qk^-Qk^ 2^ hk-k",ai • (hk"-k',a2 X hk'-k.aa) 

aba'b' 010203 kk' k" 

X Imj^ry^ G'kbai(»£p + iUJq)G]^-aia2{'i^p + ^^9)Gk'a2a' + »^g)G'k'ya3 (^ep)G'ka3a , (10) 

where we have introduced the bare Green's function, 

^ /. X \ - (ai|ua(k))(Ma(k)|a2) 

^ lEp - (ea(k) - /i) + 7Sign(ep) 

with |MQ(k)) and eQ(k) the a-th eigenfunction and eigenenergy of ■Haa'(k). l/r is a phenomenological damping rate, 
and tp = {2p + 1)ttT is a fermionic Matsubara frequency with integer p. By combining eq. ^ and (jlOp and making 
Fourier transformation back to real space, we end up with the foUowing expression 

Cxy = X! • (hj2 X hj3)W^2^(ii, 12,13), (12) 

(*l,*2,i3) 



with 



4 

T-yO 



lim— ^ (-1)-^ ^ Im|^/a;(iep, ia;g,n, 12) J(iep + ?Wg,i2,«3)-^y(«ep + iwg, --iiLjq,i3,ii) , (13) 



-P(n,i2,i3) 



where ujq is set to be ujq ~ 27rT. In eq. (|13p. the summation is taken over the 6 permutation of sites ii,i2 and 13 
with the sign of permutation (—1)^. The weighting factor Wxy{ii, 121^3) contains information about the electronic 
structure of the system through the functions = I • e^j and J: 

1 1. ^ _ ^ (14) 

J{iep;ii,i2) = ivEG'kaia2(*ep)e^^'"'i '''"S 

k 

where ai and a2 are the sublattice indices corresponding to sites ii and Z2. 



III. QUANTITATIVE CONSISTENCY WITH EXPERIMENTS 

While our approach based on the single-orbital model gives qualitative description of the Hall conductivity of 
Pr2lr207, it requires a realistic band structure to achieve a quantitative agreement with experimental data. In 
particular, the amplitude of Hall conductivity is sensitive to orbital degeneracy. In Pr2lr207, the conduction bands 
are composed of three-fold degenerate t2g orbitals of Ir 5d electrons. In our perturbative framework, orbital degeneracy 
can be incorporated by replacing the summation over sites in eq. (3) in the main text, by the summation over both 
sites and orbitals. This leads to an increase of transverse scattering channels, potentially enhancing an considerably. 

To consider the effect of orbital degeneracy, we adopt a multi-orbital tight-binding model with three t2g orbitals, 
with Slater Koster determination of transfer integrals via Ir-O-Ir hopping paths. This tight-binding model is equivalent 
to the one used in ref. [ij . We choose an electron density close to that of d^ configurations in the three-fold degenerate 
t2g orbitals expected in Pr2lr207: n = 0.85. As parameters relevant to Pr2lr207, we choose 1/t = 0.05, J = 0.073, the 
lattice parameter lOA, Jspin = 3.48K and the magnetic moment of Pr'^'*' multiplct as 3.0/iB. This choice of parameter 
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sets gives a diagonal resistivity p ~ 8.0 x lO^il ^cm ^ within the first Born approximation, of the same order as the 
experimental value ~ 2.8 x lO^f^^^cm-^d, Q. 

As a result, we obtained Hall conductivity as large as ~ 30i7~^cm~^ at high fields [H > 4.0 Tesla) for H |j [100], 
comparable to the experimental data, in which the plateau of cth is observed at H > 4.0 Tesla, taking an almost 
constant value ~ 30r2~^cm~^. Furthermore, for H || [111], the Hall conductivity has a peak around H ^ 1.0 Tesla, 
and shows sign reversal at H ^ 8.0 Tesla, again consistent with experimental data where the peak is found at 0.7 
Tesla Q, while the sign reversal takes place around H ^ 6.0 Tesla Q. 
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